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Abstract 

We show that the distribution of the square of the supremum of reflected fractional 
Brownian motion up to time a, with Hurst parameter-H greater than 1/2, is related 
to the distribution of its hitting time to level 1, using the self similarity property 
of fractional Brownian motion. It is also proven that second moment of supremum 
of reflected fractional Brownian motion up to time a is bounded above by a 2H . 
Similar relations are obtained for the supremum of fractional Brownian motion 
with Hurst parameter greater than 1/2, and its hitting time to level 1. What is 
more, we obtain an upper bound on the complementary probability distribution 
of the supremum of fractional Brownian motion and reflected fractional Brownian 
motion up to time a, using Jensen's and Markov's inequalities. A sharper bound is 
observed on the distribution of the supremum of fractional Brownian motion by the 
properties of Gamma distribution. Finally, applications of the given results to the 
financial markets are investigated and partial results are provided. 
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1 Introduction 



Stochastic Model 



The price of the stock, a volatile asset, is often described by Black-Scholes 
model which is also known as Geometric Brownian model. 
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The stochastic differential equation of this model is 



dY t = JlY t dt + aY t dW t , < t < T (1) 

where ft is a constant mean rate of return and a is a constant volatility, W 
is the standard Brownian motion and T is the time of maturity of the stock. 
When we solve this equation explicitly, we obtain 



Y t = F exp ((// - -a 2 )t + aW t ) (2) 

where Yq is the initial value of one share of stock. However, this model has 
well known deficiencies such as increments of standard Brownian motion are 
independent. This is not what we are exposed to in real life. In order to 
overcome this problem various alternative models have been suggested such 
as fractional Brownian motion, (fBm). 

Fractional Brownian motion also appears naturally in many situations. Some 
examples are, the level of water in a river as a function of time, the characters 
of solar activity as a function of time, the values of the log returns of a stock, 
the prices of electricity in a liberated electricity market, pp. In the first three 
examples fBm with Hurst parameter, H > 1/2 is used which means that 
the process is persistent. And the last example is modeled by fBm with Hurst 
parameter H < 1/2. Through out this article we are interested in the financial 
applications of fBm, specifically the model of the values of log returns of a 
stock. In other words we are interested in the fBm with Hurst parameter 
H > 1/2. The Black-Scholes model for the values of the log returns of a stock 
using fBm is given as 



Y t = Y exp (0 + fj)t + aBf), < t < T (3) 

where Yq is the initial value, r is the constant interest rate, fi is the constant 
drift and a is the constant diffusion coefficient of fBm which is denoted by 
(-Bf )t>o- 

The advantage of modeling with fBm to the other models is its capability 
of displaying the dependence between returns on different days. FBm on the 
other hand is not a semimartingale and is not a Markov process and it allows 
arbitrage, [3]. However, in order to display long range dependence and to 
have semimartingale property other models can be constructed using fBm 
as a central model, [5j. Due to the complex nature of fBm, the most useful 
and efficient classical mathematical techniques for stochastic calculus are not 
available for fBm. Therefore most of the results in the literature are given as 
bounds on the characteristics of fBm. In spite of this, fBm still possess some 



nice properties, it is a Gaussian process and it holds self similarity property. 
Most of the techniques developed for fBm are related to these properties. 
Now, there is no doubt that the investors would be interested in having some 
information on value of the supremum of the asset in order to manage the risk 
or maybe to hedge financial assets and to construct portfolios. Therefore in 
this study, mainly two new results on the supremum of fBm and an application 
of these results are introduced. 

• An identity for the distribution of square of the supremum of fBm 

• An identity for the distribution of square of the supremum of reflected fBm 

• An upper bound on the second moment of reflected fBm 

• A lower bound on the distribution of the supremum of fBm upto time 1 or 
any fixed time t. 

• A theoretical application of the above results is given, which is commonly 
used in the literature to obtain a measure of risk in finance that is called as 
"maximum drawdown." 

Fractional Brownian Motion 

Let us start with introducing fBm. FBm was first introduced within Hilbert 
space framework by Kolmogorov, [2J. It was named due to the stochastic 
integral representation in terms of standard Brownian motion, which was given 
by Mandelbrot and Van Ness, [3]. 

Let if be a constant in (0,1). FBm (B^) t > with Hurst parameter if is a 
continuous and centered Gaussian process with covariance function 

E[B?Bf] = l -(t 2H + s 2H -\t- s\ 2H ). 

For H = 1/2, fBm corresponds to a standard Brownian motion. A standard 
fBm, (B^) t >o has the following properties: 

• B* = and E[B?\ = for all t > 0. 

• B H has homogenous increments, that is B^ +s — B^ has the same law as Bf 
for s, t > 0. 

• B H is a Gaussian process and E[(Bf ) 2 } = t 2H ,t > 0, for all H e (0, 1). 

• B H has continuous trajectories. 



Stochastic integral representation 



Several stochastic integral representations have been developed for the fBm. 
For example, it is proved that the following process is a fBm with Hurst 
parameter H E (0,1), 

^iwW"'-'*" 1 *-'-'*" 1 ^' 
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+ j(t- s) H - l l 2 dW s ) (4) 
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where W t is a standard Brownian motion with Wo = considered on a prob- 
ability space (Q,J-, P). The paths of W are continuous, the increments of 
W over disjoint intervals are independent Gaussian random variables with 
zero-mean and with variance equal to the length of the interval. And V is the 
gamma function. Now let us omit the constant 1/Y(H +1/2) for simplicity 
and use the change of variable s = tu, then 



E[{B» f\ = J((t - sfi- 1 ' 2 - (s) H + - 1,2 fds 

E 

=e H J((i-u)»- 1/2 -(-u)»- 1/2 ) 2 du 



= C(H)t 



2H 



(5) 



Also 



E[{B* - B?)*\ =j((t- uf + -^ -(s- u)%-^)*d8 

E 

=e H J((t- S - u) H + - 1/2 - (-uf + - 1/2 fdu 

R 

= C(H)\t-s\ 2H . (6) 



Hence 



h r>m 
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Correlation Between Two Increments 



For H — |, the process {B^) t > Q corresponds to a standard Brownian motion, 
in which the increments are independent. For H ^ 1/2 the increments are 
not independent. By the definition of ffim, we know the covariance between 
B H (t + h)- B H {t) and B H (s + h) - B H (s) with s + h < t and t - s = nh is 

p H {n) = \h™[{n + lf H + (n- 1) 2H - 2n 2H ]. 

We observe that two increments of the form B H (t + h) — B H (t) and B H (t + 
2h) — B H (t + h) are positively correlated for H > 1/2, and they are negatively 
correlated for H < 1/2. 

Self Similarity Property 

Since the covariance function of fBm is homogenous of order 2H, fBm possess 
the self-similarity property, that is for any constant c > 0, 

{B»)t>o l = (c H B? ) t > . (8) 

Note that, by taking H — | we obtain the self similarity property of standard 
Brownian motion (W t )t>o that is 

(W at ) t > ^ (a 1 /^)^,. (9) 

Using equation it was shown in [6] that the supremum of reflected standard 
Brownian motion up to time 1 is identical in law with the the reciprocal of 
square root of hitting time of level 1. From there on, the first moment of this 
supremum was calculated using the properties of Normal distribution. 



2 Main Results 

Summary of the main results 

In this section, we prove that the second moment of reflected fBm up to time 
a is bounded above by a 2H . In our proof, as a new approach we apply the 
properties of Gamma distribution to fBm combined with the result given in 
[7] and use equation (jSJ). 

Another new result given in this section is, we provide a lower bound on the 
distribution of the supremum of fBm up to time 1 and up to fixed time t. In 
this proof, we consider fBm up to a random exponentially distributed time T 
which is independent of the process and using the self similarity property we 
obtain an upper bound on the expected value of the supremum of fBm up to 
time 1 and we combine this result with Markov's Inequality to find the lower 
bound on the distribution of the supremum. 



Let (B^) t >Q be a fBm defined on the probability space (fl, J 7 , P). And, let us 
define the reflected fBm around 0, that is denoted as (\Bf\) t >o, 



Let 7~i := inf{t > : \B^\ = 1}. In other words, T\ is the first hitting time of 
level 1 for the reflected fBm. 

Let H a := inf{t > : Bf = a} be the first hitting time of level a where 
a G R + . Let M t H := sup 0<1)<t \ B„\, that is the supremum of reflected fBm. 
And similarly let Sj* := sup 0< „ <t S^ be the supremum of fBm. 



Notation 




(10) 



In Figured!, we display sample paths of fBm and reflected ffim with Hurst 
parameters H = 0.3, H = 0.5 and H = 0.8 from left to right. These sample 
paths are generated using the Matlab code and the algorithm proposed by 
ABRY and SELLAN given in [TT] . This algorithm is a fast implementation of 
a method proposed by SELLAN [12] , which is using sum through a fractional 
wavelet basis. In this method, final paths carry both the short-term and long- 
term correlation information. For the details of simulations of fBm, one can 
also see a study by CAGLAR, given in 



Fig. 1. Fractional Brownian motion and their reflections around zero with Hurst 
parameters H=0.3,H=0.5,H=0.8 
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Theorem 2.1 For fBm with Hurst parameter H > | and a e R + 
( M f) 2 ^(-) 2 " and (S*) 2l ^(4-? H 

Ti Hi 



The second moment 



E(M?) 2 < a 2H 



Proof. 



For x > 0, by equation (jSj) we obtain 



P((-) 2H <x) = P(n>-^H) = P( sup |£f|<l) 



= P( sup \B?u\ < 1) = P( sup \B;\ <V^) = P{{Maf 

0<u<a * 1/2H 0<u<a 



Hence, it is seen that 
Now note that, 



P( sup \B*\ <^) = P( sup B u a < v^, inf B* > yfi) (12) 

0<u<a 0<u<a 0<«<a 

<P( sup B* < y/x) = P((S«) 2 < x) 

0<u<a 



iH\2 law , a ^2H 



And by following the same argument given in (TTTp it is also observed that 

(3f) S 



#1 



Hence, combining the equations (|TT|) and (|T2|) the following inequality can be 
obtained 

P((M? f <x)< P{{Sy <x) = P{{^-) 2H < x) 

Now, by applying the properties of the Gamma distribution, for x > 0, k > 
and 9 > we have 



k— 1 — H 

x e e 



cfe = (13) 



r(*) 

Therefore 

oo 

^((Mf ) 2 ) < E((S» f) = a™E((±r) = a 2H J E(e-* H ? H )dx 

1 o 

For a standard Brownian motion W, it is well known that E(e~ XTa ) = e~ a ^ 2 ~ x , 
where T a denotes the first hitting time of level a for standard Brownian motion. 
It was shown in [7] that, for ffim with Hurst parameter H > ^ the following 
inequality holds 

E{ e - XH ™) < e~ aV2X ~ (14) 

for all A, a > 0. 

Therefore, by equation (fl4l) 

oo oo 

a 2H J E(e- XT ? H )dx < a 2H J e'^dx = a 2H (15) 



As a result, E((M^) 2 ) < a 



□ 



Corollary 2.1 For fBm with Hurst parameter H > ^ and a£l 



■+ 



a 



H 



P(M a H > x) < P(S^ >x)< 



x 



Proof. 

Note that, by Jensen's inequality 

E\S») < E{{S»f) < a™ 
As a consequence, E(M?) < E(S^) < a H . 
Now, by Markov's inequality, it is seen that 

P(Mf > x) < P(S» >x)< — 



The bounds given above already are very useful information for the investors 
because it tells the investor the probability of the values of the supremum up 
to a fixed time. 

In the next theorem a closer upper bound is obtained for the distribution 
of the supremum of fBm which actually provides better information for the 
investors. 

Theorem 2.2 Consider fBm up to time a, with Hurst parameter H > |, 



Now, consider taking fBm up to time T, where T is exponentially distributed 
random variable with mean 1/A, independent of underlying fBm, B. Then 
using self similarity property of fBm, 



x 



□ 



P(S? >x)< 



V2a H 



x^/ti 



Proof. 




(16) 



and hence 




(17) 



Again by the self similarity property we have, 

E( sup B u ) = E(T 1/2 sup B s ) = E(v / T)^(S'i) 



0<»<T» 



where E{T p l 



2^ _ rifg) 

" Ap/2 • 



One can see [8], for details related to taking a random process up to a random 
time which is independent of the underlying process. 

Therefore we obtain 



And by the Markov's Inequality for any x > 0, we can write 

m > ,) < *m = 4. (19) 



Now by using the scaling property, for any fixed time a we observe 



P( sup 5^ < x) =P(a H sup <x)= P(a H S? < x) (20) 

0<u<a 0<m<1 

□ 



3 Application of the results to financial mathematics 



In finance, depending on the investor's demand, risk can be defined in many 
ways. One of the commonly used definitions of risk is known as "downfalls' 1 
or "maximum drawdown" . Maximum drawdown is the highest possible loss 
in the price of one share of the risky asset. One of the fundamental results 
related to a standard Brownian motion is Levy Theorem (see [1]). This theo- 
rem leads us to the result that " maximum drawdown" , the highest possible 
loss in the trajectories of the standard Brownian motion is identical in law 
with the supremum of the reflected Brownian motion. In other words, the 
"maximum drawdown" , before time t is defined as 

:= sup (W u - W v ) = sup ( sup W u - W v ) 

0<u<v<t 0<v<t 0<u<v 



and it satisfies 



( sup W u - W v ) l = ( sup \W V \ 

0<u<v 0<v<t 



by the Levy isomorphism. 

In the above result, we would like to replace standard Brownian motion with 
ffim, because it is obvious that ffim is a more realistic model than a standard 
Brownian motion for the risky assets. In the literature, due to the complex 
nature of fBm there are not many exact results on the distributions related 
to the trajectories. For this reason, we also expect to obtain not exact results 
but some bounds. 

In the fBm notation let, 

Dt' H := sup (B a , - B a ) = sup ( sup B a , - B a ) 

0<a'<a<t 0<a<t 0<a'<a 

be the "maximum drawdown" . 

Our goal is to find the exact distribution or alternatively find bounds for 
the distribution of "maximum drawdown" , which is directly applicable to 
finance, however we have not been able to obtain that, yet. In this study, 
we would like to present a lower bound for the distribution of the difference 
between the supremum of fBm and the value of fBm at time a, which is the first 
step of finding a lower bound for the distribution of "maximum drawdown" 
of fBm up to time t. And this following result will be extended to finding the 
distribution of "maximum drawdown" in the later studies. 

Now, let 

Y a H := flf - B; (21) 

be the difference between the supremum of fBm up to time a and the value of 
fBm at time a. 

Theorem 3.1 For y > 0, 

V2a H 



P(Y a H <y)>l 



Proof. 

For a E K+, 



P(S? >b) = P(S? > b, Bf >b)+ P(S" > b, Bf < b) 



P(Pf > b) + P(5f > 6, < 6) < ^£ (22) 



by Theorem ( 12.2ft . 
Thus we have, 

< P(S? > 6, B? < 6) < ^ - P(Pf > 6) (23) 
For b > x, x E R+, 

P{B?<x)-P{S?<b,B?<x) 
= P{S« > x, B» <x)< P(S» > b, B» <b)< ^f- - P(B* > 6) (24) 

by (|23|) and since x < b. Therefore, we have 



P(5f < 6, flf < s) > P(Pf < x) + P(Pf > 6) - ^ (25) 
Now result, we get 



P(Xa < V, B a <*)> P(Bj? <x)+ P(Pf > y + x) - {y f^ (26) 
by (12 ip , where y = b — x and 



finally 



Hence, for y > 0, 

OO OO j— 

p(y a H <y)>2 + J dp ( B « >y+x)-J v 2 







n 



dx 



1 V2a H 



OO 

+ J dP(B^ >y + x) (27) 



2 yfiiy 

And by the Fundamental Theorem of Calculus, 

1 V2a H °f 1 ~i*+v) 2 , , v 7 ^ 



H 



P(Y a H <y)>-- + -^e^^dx = 1 - (28) 



2 v^y 7 v/27r \/7r2/ 



o 

□ 



4 Conclusion 



The Black-Scholes model for the values of the log returns of a risky stock using 
fBm is a more realistic model than using standard Brownian motion, due the 
advantage of its capability of displaying the dependence of increments. The 
investors of the risky stocks or their derivatives, would naturally be interested 
in the highest value of the asset, in a certain time period. Towards that end, 
we have provided an upper bound on the second moment of the supremum 
of reflected fBm as a first result, and based on this, we have given an upper 
bound on the distribution of the supremum of reflected fBm. 

Later, we have observed a lower bound on the distribution of the supremum 
of fBm up to time 1 as well as up to fixed time a. As an application of the 
second result, we have given the joint distribution of the supremum of fBm 
up to time a and the value of fBm at time a and obtained a lower bound on 
the distribution of the difference between them. Obviously, this application 
is already very useful for the investors. Based on the given results, we will 
investigate the distribution of " maximum drawdown" because of its important 
use as a measure of risk, in finance, as future work. Our conjecture on this 
distribution is, that it is both related to the distribution of the supremum of 
the reflected fBm and the distribution of the difference between the supremum 
of fBm and the value of fBm at the maturity time. 
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